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Abstract 

The flavor violating fj,~ — e~ conversion in nuclei is studied within the minimal 
supersymmetric standard model. We focus on the R-parity violating contributions 
at tree level including the trilinear and the bilinear terms in the superpotential as 
well as in the soft supersymmetry breaking sector. The nucleon and nuclear structure 
have consistently been taken into account in the expression of the [i~ — e~ conversion 
branching ratio constructed in this framework. We have found that the contribution 
of the strange quark sea of the nucleon is comparable with that of the valence quarks. 
From the available experimental data on /x~ — e~ conversion in 48 Ti and 208 Pb and the 
expected sensitivity of the MECO experiment for Al we have extracted new stringent 
limits on the R-parity violating parameters. 



PACS: 12.60Jv, 11.30.Er, 11.30.Fs, 23.40.Bw 



KEYWORDS: Lepton flavor violation, exotic fi — e conversion in nuclei, supersymmetry, 
R-parity violation, muon capture. 



On leave of absence from the Joint Institut for Nuclear Research, Dubna, Russia 



1 Introduction 



The lepton flavor violating process of neutrinoless muon-to-electron (/i — e) conversion in a 
nucleus, represented by 

fjL- + {A,Z) — e" + {A, Z)* , (1) 

is an exotic process very sensitive to a plethora of new-physics extensions of the standard 
model(SM) In addition, experimentally it is accessible with incomparable sensitivity. 

Long time ago Marciano and Sanda Jl] has proposed it as one of the best probes to search 
for lepton flavor violation beyond the standard model. Recently, in view of the indications 
for neutrino oscillations in super-Kamiokande, solar neutrino and LSND data, new hope has 
revived among the experimentalists of nuclear and particle physics to detect other signals 
for physics beyond the SM. A prominent probe in this spirit is this exotic process (1). 

The fact that the upper limits on the branching ratio of the //~ — e~ conversion relative 
to the ordinary muon capture, 



R, e - = r( M - -> e-)/r( M - - !/„), (2) 

offer the lowest constraints compared to any purely leptonic rare process motivated a new 
fi~—e~ conversion experiment, the so called MECO experiment at Brookhaven [§,0,0, which 
got recently scientific approval and is planned to start soon. The MECO experiment is going 
to use a new very intense [l~ beam and a new detector operating at the Alternating Gradient 
Synchrotron (AGS). The basic feature of this experiment is the use of a pulsed \i~ beam to 
significantly reduce the prompt background from ir~ and e~ contaminations. For technical 
reasons the MECO target has been chosen to be the light nucleus 27 Al. Traditionally the 
(j, — e conversion process was searched by employing medium heavy (like 48 Ti and 63 Cu) 
||, HJ] or very heavy (like 208 Pb and 197 Au) ||, O, [12|] targets (for a historical review on such 



experiments see Ref. JL3|). The best upper limits on R pe - set up to the present have been 
extracted at PSI by the SINDRUM II experiments resulting in the values 

R^- < 6.1 x KT 13 for 48 Ti target [9], (3) 
R^- < 4.6 x KT 11 for 208 Pb target |TO|. (4) 

(at 90% confidence level). The experimental sensitivity of the Brookhaven experiment is 
expected to be roughly 

Rue- < 2 x KT 17 for 27 A1 target [6, [7[ (5) 

i.e. three to four orders of magnitude below the existing experimental limits. 

It is well known that the process (1) is a very good example of the interplay between 
particle and nuclear physics attracting significant efforts from both sides. The underlying 
nuclear physics of the \l — e conversion has been comprehensively studied in Refs. [|Hf|-[n>|. 
^From the particle physics point of view, processes like \i — e conversion, is forbidden in 
the SM by muon and electron quantum number conservation. Therefore it has long been 
recognized as an important probe of the flavor changing neutral currents and possible physics 
beyond the SM 
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On the particle physics side there are many mechanisms of the \i — e conversion con- 
structed in the literature (see |5|, |17], [18], |19| and references therein). All these mechanisms 
fall into two different categories: photonic and non-photonic. Mechanisms from different 
categories significantly differ from the point of view of the nucleon and nuclear structure 
calculations. This stems from the fact that they proceed at different distances and, there- 
fore, involve different details of the structure. The long-distance photonic mechanisms are 
mediated by the photon exchange between the quark and the [i — e-lepton currents. These 
mechanisms resort to the lepton-flavor non-diagonal electromagnetic vertex which is presum- 
ably induced by the non-standard model physics at the loop level. The contributions to the 
fi — e conversion via virtual photon exchange exist in all models which allow — > decay. 
The short-distance non-photonic mechanisms contain heavy particles in intermediate states 
and can be realized at the tree level, at the 1-loop level or at the level of box diagrams. 

The non-photonic mechanisms are mediated by various particles in intermediate states 
such as W, Z-bosons [0-0) Higgs bosons [0, |[, supersymmetric particles (squarks, sleptons, 
gauginos, higgsinos etc.) with and without .R-parity conservation in the vertices. 

In the supersymmetric (SUSY) extensions of the SM with conserved R-parity (i? p SUSY) 



fi — e conversion has been studied in Refs. [16|, |20J . In this case the SUSY contributions 



appear only at the loop or box level and, therefore, they are suppressed by the loop factor. 
The situation is different in the SUSY models with R-parity non-conservation (_R P SUSY). In 



this framework there exist the tree level non-photonic contributions |17|, [19| an d the 1-loop 
photonic contributions significantly enhanced by the large logarithms 



The primary purpose of this work is to offer a theoretical background for the running and 
planned fi — e conversion experiments. We consider all the possible tree level contributions 
to fi — e conversion in the framework of the minimal SUSY model with most general form of 
R-parity violation (_R P MSSM) including the trilinear ^couplings and the bilinear ^ p lepton- 
Higgs terms. We also examine some non-SUSY and i? p SUSY mechanisms previously studied 
T5| , [To] , 2~0|. We develop a formalism of calculating the fi — e conversion rate for the quark 



m 



level Lagrangian with all these terms. 

In our study we pay special attention to the effect of nucleon and nuclear structure 
dependence of the /i — e conversion branching ratio R^ e - ■ in particular, we take into 
account the contribution of the strange nucleon sea which, as we will see, gives a contribution 
comparable to the contribution of the valence quarks of a nucleon. 

Thus, we apply our formalism to the case of nuclei 48 Ti and 208 Pb by calculating numeri- 
cally the muon-nucleus overlap integral and solving the Dirac equations with modern neural 
networks techniques and using the PSI experimental data. A similar application is done for 
the 27 Al target by employing the sensitivity of the designed MECO experiment. 

Our final goal is to derive on this theoretical basis the experimental constraints on the lfi p 
Yukawa couplings, the lepton-Higgs mixing parameters and on the sneutrino VEVs. Towards 
this end we use the experimental upper limits on the branching ratio R^ e - given above and 
derive the new stringent constraints on the R-parity violating parameters. 
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2 The effective fi — e conversion Lagrangian 



It is well known that, the lepton flavors (Lj) and the total baryon (B) number are conserved 
by the standard model interactions in all orders of perturbation theory. As mentioned above, 
this is an accidental consequence of the SM field content and gauge invariance. Thus /x — e 
conversion is forbidden in the SM. In contrast to Li the individual quark flavors Bi are not 
conserved in the SM by the charged current interactions due to the presence of non-trivial 
Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix. Lj are conserved since the analogous 
mixing matrix can be rotated away by the neutrino fields redefinition. The latter is possible 
while the neutrino fields have no mass term and, therefore, are defined up to an arbitrary 
unitary rotation. As soon as the model is extended by inclusion of the right-handed neutrinos 
lepton flavor violation can occur since neutrinos may acquire a non-trivial mass matrix. 

In the supersymmetric extensions of SM Li and B conservation laws are in general vio- 
lated. As a result potentially dangerous total lepton {Tjj ) and baryon {]$ ) number violating 
processes become open. 

One may easily eliminate $ and Lf$ interactions from a SUSY model by introducing 
a discrete symmetry known as R-parity. This is a multiplicative Z 2 symmetry defined as 
R p = (—i^ 3B + L + 2S ) where S is the spin quantum number. In this framework neutrinos are 
massless. However the flavor violation in the lepton sector can occur at the 1-loop level via 
the Lj-violation in the slepton sector. Thus, — e~ conversion is allowed in the SUSY 
models with R-parity conservation. 

There is no as yet convincing theoretical motivation for R-symmetry of the low energy 
Lagrangian and, therefore, SUSY models with (i? p SUSY) and without (Ifl p SUSY) R-parity 
conservation are "a priori" both plausible. 

Despite the above mentioned problems with L 1 , Lfi interactions Jji v SUSY looks rather 
attractive, since it may offer a clue to the solution of some long standing problems of particle 
physics, such as neutrino mass problem. In the 1/ip SUSY framework neutrinos acquire 
Majorana masses at the weak-scale via mixing with the gauginos and higgsinos as well as 



via Tjj loop effects [21]. Furthermore, lfl p SUSY models admit non-trivial contributions to 
the lepton flavor violating processes. During the last few years the Tjl v SUSY models have 
been extensively studied in the literature (for a recent review see Refs. |J22| ) . 

We analyze possible mechanisms for process (1) existing at the tree level in the minimal 
Ijip SUSY model with a most general form of R-parity violation. 

A most general gauge invariant form of the R-parity violating part of the superpotential 
at the level of renormalizable operators reads 

W % = \ ]k L % L 3 El + X'^QjDt + HjLjH 2 + ~K jk U*Dp c k , (6) 

where L, Q stand for lepton and quark doublet left-handed superfields while E c , U c , D c for 
lepton and up, down quark singlet superfields; Hi and H 2 are the Higgs doublet superfields 
with a weak hypercharge Y = —1, +1, respectively. Summation over the generations is 
implied. The coupling constants A (A") are antisymmetric in the first (last) two indices. The 
bar sign in A', A" denotes that all the definitions are given in the gauge basis for the quark 
fields. Later on we will change to the mass basis and drop the bars. Henceforth we set A" = 
which are irrelevant for our consideration. This condition ensures the proton stability and 
can be guaranteed by special discreet symmetries other than R-parity. 
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The R-parity non-conservation brings into the SUSY phenomenology the lepton number 
[Jjj ) and lepton flavor ) violation originating from the two different sources. One is 
given by the Ijj trilinear couplings in the superpotential W$ p o{ Eq. (||). Another is related 
to the bilinear terms in W^ p and in soft SUSY breaking sector. Presence of these bilinear 
terms leads to the terms linear in the sneutrino fields i>i in the scalar potential. The linear 
terms drive these fields to non-zero vacuum expectation values (z/j) ^ at the minimum of 
the scalar potential. At this ground state the MSSM vertices Zv v and Wev produce the 
gaugino- lepton mixing mass terms Zv(y), We(y) (with W, Z being wino and zino fields). 
These terms taken along with the lepton-higgsino ^LiHi mixing from the superpotential 
of Eq. (^) form 7x7 neutral fermion and 5x5 charged fermion mass matrices. For the 
considered case of [i — e conversion the only charged fermion mixing is essential. The charged 
fermion mass term takes the form 

£$L = -^( T )-M±^' (+) - H.c. (7) 

in the basis of two component Weyl spinors corresponding to the weak eigenstate fields 

= (e~ L , fT L , T~ L , -zA_, HD, = (e+ , //+ , r+ , H+). (8) 

Here \± are the SU 2 l gauginos while the higgsinos are denoted as Hf 2 . These fields are 
related to the mass eigenstate fields *f?(±) by the rotation 

* (±)i = a± *; ±)i , (9) 

The unitary mixing matrices A^ diagonalize the char gino- charged lepton mass matrix as 

(A-)*M± (A+) T = Diag{mf\m x ±}, (10) 

where mf and m x ± are the physical charged lepton and chargino masses. In the present 
paper we use the notations of Refs. p3| |4 |. 



Rotating the MSSM Lagrangian to the mass eigenstate basis according to Eq. (||) one 
obtains the new lepton number and lepton flavor violating interactions in addition to those 
which are present in the superpotential Eq. (^|). Note that the mixing between the charged 
leptons (e + L , n + L ) and the chargino components (— H2), described by the off diagonal 
blocks of the A + , is proportional to the small factor ~ m e ^/ Msusy [0, and is, therefore, 
neglected in our analysis. 

Let us write down the MSSM terms generating by the rotation (|9[) the new lepton flavor 
violating interactions relevant for the — e~ conversion. In the two component form they 



can be written as 25 



C «SSM = ^h~ Z ^' iAi ^ + ^2A-W2, (11) 

where A^j = (1 — 2 sin 2 9w)8ij + 5^5^. Rotating this equation to the mass eigenstate basis 
we write down in the four- component Dirac notation 

C mssm = 77—2 azZ^ej^n + g 2 (i ■ u k P R e\d Lk . (12) 

L COS u\y 
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Here = (e, /i, r), P L R = (1 7s)/2. The lepton flavor violating parameters in this formula 
are given by 

a z = A u A-* 24 « £ u &, Ci = A~; 4 « (13) 



The approximate expressions for these parameters were found by the method of Ref. [23 | 
expanding the mixing matrix A~ in the small matrix parameter 

" y/2 {M211 - sin2/3M^) ' 1 j 

where /i, M 2 are the ordinary MSSM mass parameters from the superpotential term itHiH 2 
and from the SU 2 gaugino soft mass term M 2 W k W k . The MSSM mixing angle is defined as 
tan/? = (H 2 ) / (Hi). The other lepton flavor violating interaction contributing to the /i~ — e~ 
conversion come from the superpotential (|6|). 

The leading diagrams describing possible tree level Iflp MSSM contributions to the \l — e 
conversion are presented in Fig. 1. The vertex operators encountered in these diagrams are 



£ Me - = 2 Aai " Li eP L fi + Kjj "lA P l d 3 ~ x '*jk V nj («L e * P r d * + d ^ u ^ P r e i) + 
+92 Q V nk ■ u n P R e c t d Lk + l-^-azZ^P.ii - -^—Z^q lt ,(e L P L + e R P R )q, (15) 

Z COS v\v cos "W 

where the first three terms originate from the superpotential the fourth and fifth terms 
correspond to the chargino-charged lepton mixing terms in Eq. fll2|) and the last one is the 
ordinary SM neutral current interaction. 

In Eq. fll5|) the SM neutral current parameters are defined as usual 

e L = T 3 - sin 2 9 W Q , e R = — sin 2 6 W Q 

with T 3 and Q being the 3rd component of the weak isospin and the electric charge. 

The Lagrangian ([TBI) is given in the quark mass eigenstate basis which is related to the 
flavor basis q° through 

Ql,r = Vl R ■ q°. 
For convenience we introduced the new couplings 



Kjk ~ Kmn iyf) jm ( V R) kn - 



The CKM matrix is defined in the standard way as V = V£V?. 

Integrating out the heavy fields from the diagrams in Fig. 1 and carrying out Fierz 
reshuffling we obtain the 4-fermion effective Lagrangian which describes the fi — e conversion 
at the quark level in the first order of perturbation theory. It takes the form 



G 



^J2 J ^ Y' L "uL(i) ^ 'IR u uR{i) ^ 'IL °dL(i) ^ 'IR dR(i)\ 1 fa 



G 



F 



(16) 
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Figure 1: Leading l/i p MSSM diagrams contributing to y, — e conversion at the tree level, (i) 
The upper diagrams originate from the trilinear A, A' terms in the superpotential Eq. (|6|). 
(ii) The lower diagrams originate from the chargino-charged lepton mixing schematically 
denoted by crosses (X) on the lepton lines. 



The index i denotes generation so that = u,c,t and dj = d,s,b. The coefficients r] 
accumulate dependence on the Ifip SUSY parameters in the form 



1 X' X'* A M 2 

V? = ~ E V u^m +2(1-1 sin 2 M a z + 4^ CiQ^\V m \ 2 , 



v d R = 4 E r*t*™ V ™ Vnl + \ sin2 6w az > (17) 

2 8 
rjf = -2(1- - sin 2 6 w )a Zl r% = --sin 2 6 W a z 



^ G F ml {n) ' ^ G Fml (n) 



Here ro g ( n ),m v ( n ) are the squark and sneutrino masses. In Eq. (|16|) we introduced the color 
singlet currents 

where % = (u u di). 

Since in the next sections we report the new results for the nuclear matrix elements of 
48 Ti, 27 A1 and 208 Pb we are going also to update yT — > constraints on the effective lepton 
flavor violating parameters corresponding to certain non-SUSY and i? p SUSY mechanisms 
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Figure 2: Photonic and non-photonic mechanisms of the \x~ — e~ conversion within some 
extensions of the standard model: (a-c) the SM with massive neutrinos and (a-c) R-parity 
conserving supersymmetric extensions of the SM. 



shown in Fig. 2. These mechanisms were previously studied in Refs. [T6|, |20| . Let us shortly 
summarize these mechanisms for completeness. The long-distance photonic mechanisms 
mediated by the photon exchange between the quark and the [i — e-lepton currents is realized 
at the 1-loop level as the v — W loop [Fig. 2(a)] with the massive neutrinos V{ and the loop 
with the supersymmetric particles such as the neutralino(chargino)-slepton(sneutrino) [Fig. 
2(d)]. In the R-parity violating SUSY models there are also lepton-slepton and quark-squark 
loops generated by the superpotential couplings XLLE C and X'LQD C respectively |T3. The 
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short-distance non-photonic mechanisms in Fig. 2 contain heavy particles in intermediate 
states and is realized at the 1-loop level [Fig. 2(a,b,d,e)] or at the level of box diagrams [Fig. 
2(c)]. The 1-loop diagrams of the non-photonic mechanisms include the diagrams similar to 
those for the photonic mechanisms but with the Z-boson instead of the photon [Fig. 2(a,d)] 
as well as additional Z-boson couplings to the neutrinos and neutralinos [Fig. 2(b,e)]. The 
box diagrams are constructed of the W-bosons and massive neutrinos [Fig. 2(c)] as well 
as similar boxes with neutralinos and sleptons or charginos and sneutrinos. The branching 
ratio formula for these mechanisms is given in the next section. 

3 Nucleon and nuclear Structure dependence of the 
/i — e conversion rates. 

One of the main goals of this paper is the calculation of the /i — e conversion rate using 
realistic form factors of the participating nucleus (A,Z). This can be achieved by applying 
the conventional approach based on the well known non-relativistic impulse approximation, 
i.e. treating the nucleus as a system of free nucleons f|. To follow this approach as a first 
step one has to reformulate the \l — e conversion effective Lagrangian (^) specified at the 
quark level in terms of the nucleon degrees of freedom. 

The transformation of the quark level effective Lagrangian, £e//> to the effective La- 
grangian at the nucleon level, jC^ ff , is usually done by utilizing the on-mass-shell matching 



condition \2t 



(y F \£l ff \yj) « OMC^I*,), (19) 

where and (^f\ are the initial and final nucleon states. In order to solve this equation 
we use various relations for the matrix elements of the quark operators between the nucleon 
states 

(N\q T K q\N) = G { £ N) V N T K V N , (20) 

with q = {u,d,s}, N = {p,n} and K = {V, A, S, P}, T K = {7^,7^75,1,75}. Since the 
maximum momentum transfer q 2 in fi — e conversion, i.e. |q| ~ m^/c with = IQb.QMeV 
the muon mass, is much smaller than the typical scale of nucleon structure we can safely 
neglect the q 2 -dependence of the nucleon form factors G% . For the same reason we 
drop the weak magnetism and the induced pseudoscalar terms proportional to the small 
momentum transfer. 

The isospin symmetry requires that 

r*(u,n) _ n (d,p) _ n d ri( d ' n ) — r<( u 'P) — r<u r<( s > n ) — n( s 'P) — n** fot\ 

with K — V,A, S, P. Conservation of vector current postulates the vector charge to be equal 
to the quark number of the nucleon. This allows fixing of the vector nucleon constants 

G«=2, G d v = l, G s v = 0. (22) 

The axial- vector form factors Ga can be extracted from the experimental data on polarized 



nucleon structure functions [27, 28] combined with the data on hyperon semileptonic decays 
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29| . The result is 



GPi « 0.78, G\ 



-0.47, 



-0.19. 



(23) 



The scalar form factors are extracted from the baryon octet B mass spectrum Mb in com- 
bination with the data on the pion-nucleon sigma term 

a nN = (1/2) (m u + m d )(p\uu + dd\p). (24) 

Towards this end we follow the QCD picture of the baryon masses which is based on the 



relation 30, 31 



(B\0£\B) = M B BB (25) 
and on the well known representation for the trace of the energy-momentum tensor |Kj 

0£ = m u uu + m d dd + m s ss - (ba s /87r)G a ^G^ u , (26) 

where G°L, is the gluon field strength, ot s is the QCD coupling constant and b is the reduced 
Gell-Mann-Low function with the heavy quark contribution subtracted. Using these relations 
in combination with SU(3) relations |31| for the matrix elements (B\8£\B) as well as the 
experimental data on Mg and a^N we derive 



G u s « 5.1, 



G% « 4.3, G'g ss 2.5. 



(27) 



The nucleon matrix elements of the pseudoscalar quark currents can be related to the 
divergence of the baryon octet axial vector currents Utilizing this fact we find the 

pseudoscalar form factors 



Gp « 



103, G^ ~ 



yd 



100, 



G 8 P P3 3.3. 



(28) 



Note that the strange quarks of the nucleon sea significantly contribute to the nucleon 
form factors Ga, Gp and Gs- This result dramatically differs from the naive quark model 
and the MIT bag model where G S A s P = 0. The contribution of the strange nucleon sea will 
allow us to extract additional constraints on the second generation l/i p parameters. 

Now we can solve Eq. flT9| ) and write the effective fi — e conversion Lagrangian at the 
nucleon level as 



Gp 
V2 



ejfi(l - 75) )U • + e/x • J + + e7 5 /i • J~ 



where we have introduced the nucleon currents 



NY (4 0) + < ; r 3 ) + « + <r 3 ) 75 



.,(3) 



(0) 



.,(3), 



N 



a ±s + a ±'s r 3) + («±p + «±pT3)7 



.,(3) 



(0) 



.,(3) 



N, 



(29) 



(30) 
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for nucleon isospin doublet N T = (p,n). The coefficients in Eq. ( p0|) are defined as 

af = i(G^ + G$)(t# + ijf + i# + i£), 
4 3) = l(G u v - G d v )( V f +r,f - V f -rftl 
«i 0) = kG^ + G^rjf-vf + rjf-v^ + lG^vf-vf), 



4 



«? } = liGl-GiXrjf-vt-rjf + rji 1 ] 



«S = ^(G u s + G d s M 1 ±r ] d R 1 ) + ±G s s (v? ±fjf), (31) 
«S = -^(GS-G5)(^f±^). 

"I'p = ^(Cp - Gp)(f}f- TVr), 



Starting from the Lagrangian ( p9[ ) it is straightforward to deduce the formula for the total 
fi — e conversion rate. In the present paper we focus on the coherent process, i.e. ground 
state to ground state transitions. This is a dominant channel of fi — e conversion which, in 
most of the experimentally interesting nuclei, exhausts more than 90% of the total — > e~ 
branching ratio [JT1, [Tj|] . To the leading order of the non-relativistic expansion the coherent 
— e conversion rate takes the form 

r?- = %^ Q(M P + M n f , (32) 



where 



r> ol (°) i (3) i|2 | I (0) , (3) ,,2 , | (0) . (3) 

Q = 2\a v ' + ay (p\ + \a\' s + a\' s <p\ + \al' s + a_ s 

„(°) _L „(3) ^U(°) _L , ^,( 3 ) _L «(3) 



+2 Re{(c^ + a$> <j>)[a^ 8 + a™ 8 + (a$ + a™ s ) 0]} . (33) 

The transition matrix elements M. p<n in Eq. (^) depend on the final nuclear state pop- 
ulated during the /i — e conversion. We should stress that, after computing the nuclear 
matrix elements M. p>n the data provide constraints on the quantity Q of Eq. (|33|). For 
the ground state to ground state transitions in spherically symmetric nuclei the following 
integral representation is valid 

M P:n = 4tt J J (p e r)$ P (r)p P:n (r)r 2 dr (34) 

where jo(x) the zero order spherical Bessel function and p p , n the proton (p), neutron (n) 
nuclear density normalized to the atomic number Z and neutron number N of the partic- 
ipating nucleus, respectively. The space dependent part of the muon wave function ^ is 
a spherically symmetric function which in our calculations (see Sect. 4) was obtained by 
solving numerically the Shrondinger and Dirac equations with the Coulomb potential. 
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In denning Q, Eq. (|33|), we introduced the ratio 



<j>=(M p -M n )/(M p + M n )a(A-2Z)/A, (35) 

where A and Z are the atomic weight and the total charge of the nucleus. The quantity Q 
depends weakly on the nuclear parameters through the factor 0. In fact, the terms depending 
on (p are small since <fi < 1 (see Table 1) and Gg, G d s as well as G v , G v have the same sign. In 
practice the nuclear dependence of Q can be neglected and thus, Q coincides with the value 
of 2p where p is defined below. It can be considered as a universal effective 1/ip parameter 
measuring the 1/ip SUSY contribution to the p — e conversion. It also represents a suitable 
characteristic which allows comparison of p — e conversion experiments on different targets 
treating the corresponding upper bounds on Q as their sensitivities to the Ifi p SUSY signal. 

For completeness, in Sect. 4 the limits for some non-SUSY [Tj| ^] as well as i? p SUSY 
contributions to the p~ — e~ nuclear conversion (see in Fig. 2) are updated. The correspond- 



ing expression for i?» e - is written as |TJ 



R^e- = PI, ( 36 ) 



where p is nearly independent of nuclear physics JT5| and contains the lepton flavor violating 



parameters corresponding to the contributions in Fig. 2. Thus, e.g. for photon-exchange 
mode p is given by 

I a \2 I/a/1 + /so| 2 + I/El + /moP / q „n 

P = (47m) jG^tf (37) 

where the four electromagnetic form factors Jeo, Iei, /mo, /mi are parametrized in a specific 
elementary model flT5| . 

The factor 7 (A, Z) in Eq. (|36|) accumulates about all the nuclear structure dependence 
of the branching ratio R^- . Assuming that the total rate of the ordinary muon capture is 
given by the Goulard-Primakoff function, fop, the nuclear structure factor j(A, Z) takes the 
form 

1[A,Z) = 7 = 38 

m i G 2 Zf GP (A,Z) 

where G 2 ~ 6. Thus, a non-trivial nuclear structure dependence of the pT — > e~ conversion 
branching ratio R^ e - is mainly concentrated in the nuclear matrix elements M 2 [[|. In the 
proton-neutron representation one can write down 

M 2 = [M P + Q M n } 2 (39) 

where Q takes the values of Eq. (32) of Ref. [15a] and M p>n are given by writing the matrix 
elements of Eq. ( |3l| ) in terms of an effective muon wavefunction as 

M p , n = (^}M p , n (40) 

In our present approach the role of p in Eq. (|36|) is played by Q, since Q = 2p, and the 
corresponding 7 function defined in Eq. ([38]) for R-parity violating interactions, 7^ p , is 
obtained from Eq. (|38"D by putting Q = 1 in Eq. (|39|). The separation of nuclear physics 
from the elementary particle parameters is not complete but we have seen that <fi is quite 
small. In any case we present in Table 1 the values of for the various nuclear systems. 
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4 Results and Discussion. 



The pure nuclear physics calculations needed for the p, — e conversion studies refer to the 
integrals of Eq. (|34T). The results of M. p and M. n for the currently interesting nuclei Al, 
Ti and Pb are shown in Table 2. They have been calculated using proton densities p p 
from the electron scattering data [p2| and neutron densities p n from the analysis of pionic 



atom data [33]. We employed an analytic form for the muon wave function ^(V) obtained 
by solving the Schrodinger equation using the Coulomb potential produced by the charge 
densities discussed above. This way the finite size of a nucleus was taken into consideration. 
Moreover, we included vacuum polarization corrections as in Ref. |~4]]. In solving the 
Schrodinger equation we have used modern neural networks techniques which give the 
wave function ^(r) in the analytic form of a sum over sigmoid functions. Thus, in Eq. 
(|3"4|) only a simple numerical integration is finally required. To estimate the influence of the 
non-relativistic approximation on the muon wave function ^(r), we have also determined 
it by solving the Dirac equation. The results for A4 p , n do not significantly differ from those 
of the Schrodinger picture. In Table 2 we also show the muon binding energy e& (obtained 
as output of the Dirac and Schrodinger solution) and the experimental values for the total 
rate of the ordinary muon capture taken from Ref. ||35|| . 

Using the values of A4 P , M. n for a set of nuclei throughout the periodic table we can 
estimate the nuclear structure dependence of the p — e conversion branching ratio, i.e. the 
function ^y(A, Z) in Eq. (j3~6|). The results are quoted in Table 3. For comparison in this 
table we also list the results of corresponding to the Ifi p SUSY mechanisms studied in 
the present paper and calculated as we stated in Sect. 3. These results can be exploited for 
setting constraints on the quantities p and Q corresponding to specific models predicting the 
p — e process. 

In Table 4 we quote the upper bounds for the quantities p and Q corresponding to the 
mechanisms shown in Figs. 1,2. These bounds were derived from the recent experimental 
upper bounds on the branching ratio -R Me - f° r Ti and Pb targets given in (||) and (f|) and 
from the expected experimental sensitivity (|5|) of the Brookhaven MECO experiment. The 
limits on p and Q quoted in Table 4 for 27 Al are improvements by about four orders of 
magnitude over the existing ones. 

We should stress that limits on the quantities p of Eq. (^) and Q of Eq. (|33|), are 
the only constraints imposed by the p — e conversion on the underlying elementary particle 
physics. One can extract upper limits on the individual lepton flavor violation parameters 
(ftp couplings, effective scalar and vector couplings, neutrino masses etc. [[!], |5], [T^, |TJ|) under 



certain assumptions like the commonly assumed dominance of only one component of the 
p — e conversion Lagrangian which is equivalent to constrain one parameter or product of 
the parameters at a time. Using the upper limits for Q given in Table 3 we can derive 
under the above assumptions the constraints on ctjp of Eq. (p3|) and the products of various 
ft p parameters. Thus, the bounds obtained for the scalar current couplings a± s in the im- 
parity violating Lagrangian for the 27 A1 target [| are |e*S| < 7 X 1(T 10 . The limit for a^l 
obtained with the data for the Ti target || is \a±s\ < 1.1 X 1CU 7 , i.e. more than two orders 
of magnitude weaker than the limit of 27 Al. 

With these limits it is straightforward to derive constraints on the parameters of the 
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initial Lagrangian flT5|). In Tables 5,6 we list the upper bounds on the products of the 
trilinear ft p couplings which we obtained from the experimental limit on n — e conversion 
in 48 Ti and from the expected experimental sensitivity of MECO detector using 27 Al as a 
target material. The corresponding constraints for 208 Pb are significantly weaker and not 
presented here. In Tables 5,6 the quantity B denotes a scaling factor defined as 

Bn = (i^/6.1 x 1(T 13 ) 1/2 , B Al = (R e * p /2.0 x KT 17 ) 1 / 2 , (41) 

which can be used for reconstructing the limits on the listed products of the Tjl p parameters 
corresponding to the other experimental upper limits on the branching ratio R e ^P. In the 
2nd column of Tables 5,6 we present previous limits existing in the literature and taken 



from [18] , |22| . As seen from Tables 5,6 the /x — e-conversion limits on the products A' A' and 
AA', except only A / 232 A / 1 32 and A^A'^, are significantly more stringent than those previously 
known in the literature. The two products A'A' and AA' are less stringently constrained by 
the present experimental data on 48 Ti within the tree level non-photonic mechanism. Note 
that the corresponding previous constraints on these products (2nd column of Table 5) were 
obtained from the photonic 1-loop mechanism |jTB[ of yT — e~ conversion which are better 
than existing in the literature non-/i~ — > e~ constraints on A'A' and AA'. The products in 



Table 6 are not constrained by this mechanism [113]. However it constrains the products AA 



not constrained by the tree level ft p SUSY mechanism. At present the \x~~ — > e~ constraints 
on AA within the 1-loop ft p SUSY mechanism are weaker than those derived from the other 
processes. As we have mentioned at the beginning, significant improvement on these yT — > e~ 
limits is expected from the ongoing experiments at PSI || and even better from the MECO 
experiment at Brookhaven || . This would make the \x — e conversion constraints on all the 
products of the ft p trilinear coupling attainable in the /i — e conversion at the tree and at 
the 1-loop levels better than those from the other processes in all the cases. 

Note that the last four limits for A'A in Table 6 originate from the contribution of the 
strange nucleon sea. These limits are comparable to the other \l~ —* e~ constraints derived 
from the valence quarks contributions. 

Finally we put constraints on the products of the bilinear ft p parameters evaluated from 
the fx — e conversion in 48 Ti as 

(PiH^hfim, (*iK fa) Hi < (80MeV) 2 (j^^B, (42) 
where B is a scaling factor defined in Eq. (4"T). These constraints are weaker than those 



derived in Ref. [36 from the Super-Kamiokande atmospheric neutrino data. 



The limits in Tables 5,6 and in Eq. ( |42|) were extracted under the following assumptions. 



We assumed that all the scalar masses in Eq. (|ITD are equal rh U L(n) ~ rhdLMn) ~ m u{n) ~ m 



and also that there is no significant compensation between the different terms contributing 
to the ratio Rue-- Note that the last assumption is, in practice, equivalent to the other well 
known assumptions about the dominance of only one parameter or product of the parameters 
at a time. These assumptions are widely used in the literature for derivation of constraints 
on the R-parity violating parameters. 

Instead of considering only one specific combination of A and A' as dominant, one may 
attempt to consider all combinations of the R-parity violating couplings. To this end one 



13 



needs constraints on the relative magnitudes of the R-parity violating couplings. This can 
be accomplished in a fashion analogous to the theory of textures for the Yukawa couplings. 
This theory, which expresses the entries of the down fermion mass matrix in powers of a 
parameter e = 0.23, has been successful in describing the charged fermion mass matrices 
(for review see fl37] and references therein). This can be extended to the R-parity violation 



couplings themselves [38, 39]. In fact using solution B of Ref. [39| we obtain 



X' - X' - X' — X' 

/ >112 — — ^919 — 



v 121 



212 



A 



121 



3 A' 



(43) 



This way, p — e conversion can be used to constrain the overall scale of the R-parity violating 



interaction using phonomenologically acceptable textures [39] 



5 Summary and Conclusions. 

The transition matrix elements of the flavor violating p~ — e~ conversion are of notable 
importance in computing accurately the corresponding rates for each accessible channel of 
this exotic process. Such calculations provide useful nuclear-physics inputs for the expected 
new data from the PSI and MECO experiments to put more severe bounds on the muon- 
number violating parameters determining the effective currents in various models that predict 
the exotic p~ — > e~ process. 

We developed a systematic approach which allows one to calculate the p~~ — e~ conversion 
rate in terms of the quark level Lagrangian of any elementary model taking into account the 
effect of the nucleon and nuclear structure. Our conversion rate formula (^) is valid for 
the interactions with the (axial-)vector and (pseudo-) scalar quark and lepton currents, as 
shown in Eq. (|29|) . In the previous p~ — > e~ calculations found in the literature only the 
(axial-) vector currents were considered. 

In the case of the R-parity violating interactions discussed here we have investigated 
all the possible tree level contributions to the p — e conversion in nuclei. We found new 
important contribution to pT — conversion originating from the strange quark sea in the 
nucleon which is comparable with the usual contribution of the valence u, d quarks. 

We introduced the quantities p and Q defined in Eqs. (|37| ) and ( |33| ) which can be 
associated with theoretical sensitivities of a p~ — e~ conversion experiment to the charged 
lepton flavor violating interactions discussed in the present paper. These quantities are 
independent of a target material and, therefore, might be helpful for comparison of searching 
potentials of different p — e conversion experiments and for planning future experiments. 

^From the existing data on R pe - in 48 Ti and 208 Pb and the expected sensitivity of the 
designed MECO experiment |J we obtained stringent upper limits on the quantities p and 
Q. Then we extracted the limits on the products of the trilinear R p parameters of the type 
AA', A'A' which are significantly more severe than those existing in the literature. 

Let us conclude with the following important remark. As was observed in Refs. ||i~8|, |i~9"f, 



if the ongoing experiments at PSI JT(J and Brookhaven [y] will have reached the quoted 
sensitivities in the branching ratio R^ e - then the p~ — *> e~ constraints on all the products of 
the $p parameters AA, A'A, A'A' measurable in p — e conversion will become more stringent 
than those from any other processes. This is especially important since no comparable 
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improvements of the other experiments probing these couplings is expected in the near 
future. 
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0.000 


24. 


12. 


0.014 


0.000 


27. 


13. 


0.000 


-0.037 


32. 


16. 


0.023 


0.000 


40. 


20. 


0.037 


0.000 


44. 


20. 


-0.063 


-0.091 


48. 


22. 


-0.083 


-0.083 


63. 


29. 


-0.056 


-0.079 


90. 


40. 


-0.054 


-0.111 


112. 


48. 


-0.108 


-0.143 


208. 


82. 


-0.152 


-0.212 


238. 


92. 


-0.175 


-0.227 



Table 1: The variation of the quantity <fi of Eq. ( p3|) through the periodic table. For 
comparison its approximate expression m (A — 2Z)/A is also shown. 



Nucleus IpeK/m- 1 ) e b (MeV) (xlO 6 s' 1 ) A^ p (/m~ 3 / 2 ) M n (fm- 3 / 2 ) 
0.531 -0.470 0.71 0.047 0.045 

0.529 -1.264 2.60 0.104 0.127 

0.482 -10.516 13.45 0.414 0.566 

Table 2: Transition matrix elements (muon- nucleus overlap integrals M. v<n of Eq. (|34])) 
evaluated by using the exact muon wave function obtained via neural networks techniques. 
Other useful quantities (see text) are also included. 
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208 
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Mechanism 



27 Al A8 Tl 208p 6 



Photonic (7) 


4.3 


9.4 


17.3 


W-boson exchange (7) 


34.2 


25.3 


49.2 


SUSY s-leptons (7) 


11.3 


25.6 


49.5 


SUSY Z-exchange (7) 


27.6 


110.6 


236.2 


SUSY 


16.5 


46.4 


96.9 



Table 3: The nuclear physics part of the branching ratio R^, i.e. values of the function 
7$ p > 7 of Eq. (|3^) for the nuclear targets 27 Al, 48 Ti and 208 Pb. 



Mechanism 21 Al A8 Tl 208 Pb 



Photonic 


P 


< 


4.6 


xlO" 


-18 


P 


< 


8.2 


xlO' 


-14 


P 


< 3.2 


xlO" 


-12 


W-boson exchange 


P 


< 


5.8 


xl(T 


-19 


P 


< 


3.0 


xlO' 


-14 


P 


< 1.1 


xlO' 


-12 


SUSY sleptons 


P 


< 


1.8 


xl(T 


-18 


P 


< 


3.0 


xlO' 


-14 


P 


< 1.1 


xlO" 


-12 


SUSY Z-exchange 


P 


< 


7.3 


xlO" 


-19 


P 


< 


0.7 


xlO' 


-14 


P 


< 0.2 


xlO' 


-12 


R p SUSY 


Q < 2.6 


•10" 


-19 


Q < 0.7 -10" 


-14 


Q < 1.1 


•10" 


-13 



Table 4: Experimental limits on the elementary particle sector of the exotic fi — e conversion. 
Quantities p and Q are defined in of Eqs. (|37| ) and Eq. (|33|), respectively. The limits are 
extracted from the recent experimental data for the nuclear targets 48 Ti and 208 Pb |10| , p~T|| 
and from the expected sensitivity of the MECO experiment for the 27 Al target 0[see Eqs. 

(!)-(!)]• 
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Table 5: Upper bounds on A^- fc A; mn derived in the tree-level ljl v SUSY mechanism from the 
SINDRUM II data on the /i~ — e~ conversion in 48 Ti [Eq. (0)] and from the expected 
sensitivity [Eq. (§)] of the MECO experiment at BNL with 27 A1. For comparison we present 
previous bounds derived in the 1-loop Tjl v SUSY mechanism JTB[. In brackets [ ] we show 



bounds from the other processes (for references see Refs. [22, [18[]) when they are more 



stringent then the 1-loop /i — e conversion bounds. The limits are given for the scalar 
superpartner masses fh = 100 GeV. The scaling factor B is defined in Eq. (|41|). 
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2.6 ■ 10~ n 


|A 3U A321 


6.0 • 10~ 3 


4.1 • 10~ 9 


2.6 • 10" 11 


A' 2 22 A212 


9.0 • 10- 3 


7.7- 10~ 9 


5.0- 10- 11 


A322 A312 


1.2 • 10~ 2 


7.7- 10- 9 


5.0- 10- 11 


A' 12 2 A 121 


1.0 • 10- 3 


7.7- 10~ 9 


5.0- 10- 11 


A322 A321 


1.2 ■ 10~ 2 


7.7- 10- 9 


5.0- 10- 11 



Table 6: The same as in Table 5 but for the upper bounds on A^- fc 
are taken from the other processes (for references see Refs. |^2|. 1 



\imn- Previous limits 
The 1-loop $ p SUSY 



mechanism [Tj|] does not constraint these products of the lfi p couplings. 
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